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A.   Problem  Statement 

Bearings-only  data  for  a  single  target  from  two  sensors 
which  may  be  moving  or  stationary  are  available  at  two  distinct 
times.   The  following  quantities  are  required:  an  estimate  of 
course,  speed  and  position  of  the  target  at  the  latest  time; 
an  estimate  of  a  future  position  of  the  target  and/or  an  estimate 
of  a  point  on  the  track  of  the  target  with  a  specified  lead  distance 
at  a  future  time.    The  relative  positions  of  the  two  sensors 
are  assumed  to  be  known  at  the  time  of  each  target  bearing 
determination . 


B.   Operational  Analysis 

Two  simultaneous  bearings  from  two  sensors  at  two 
distinct  times  and  with  known  relative  positions  are  used  to 
estimate  the  course  and  speed  of  a  target.   The  HP-67/97  program 
presented  here  was  designed  so  that  the  data  corresponding  to 
the  earliest  time  point  is  purged  if  data  corresponding  to  a 
third  time  point  is  introduced.   The  relative  position  of  the 
sensors  may  be  updated  when  required.   Thus  the  estimated  target 
position,  course  and  speed  are  continually  updated  as  new 
information  becomes  available.   No  course  smoothing  is  performed. 


C.   Computational  Algorithm 

1.  Enter  the  course  \p        and  speed  V   of  the  primary  sensor 

V 

2 .  Enter  the  bearing  <$>      and  range   p   of  the  secondary  sensor 
Sp   from  the  primary  sensor   S,   at  the  time  of  the  latest 
bearing  observation. 

3.  Enter  the  time   t, ,  the  bearing  of  the  target  from   S, ,  and 
the  bearing  of  the  target  from   S^.   Output  the  target  range 
from   S,  . 

4.  Repeat  Step  3  or  Steps  2  and  3  for  a  second  time   t«  >  t, . 

5.  Compute  and  output: 

a.  The  estimated  course  and  speed  of  the  target. 

b.  The  bearing  and  range  (n.mi.)  of  the  target  from   S, 
at  time  t~ . 

c.  The  bearing  and  range  (n.mi.)  of  the  target  from   S_ 
at  time   t„. 

6.  If  required,  enter  a  time  tg    >  t«   at  which  a  lead  distance 
I       (n.mi.)  is  required.   Then  compute  and  output  the 
target's  predicted  bearing  and  range  from  both   S,   and   S?. 

7.  Repeat  from  Steps  1,  2,    3  or  4  as  required. 


D.   HP-67/97  Calculator  Program 
1.   User  Instructions 


Instruction 


Input 


Key(s) 


Output 


Enter  program  card  (both  sides) . 

Set  for  HP-67  output  (default)  or 
set  for  HP-97  output. 

Enter  course  and  speed  of  the 
primary  sensor   S, . 

Enter  bearing  and  range  of   S 
from   S..   at  datum. 


Enter  time   t.   of  datum, 
i 

bearing  from  S  ,  and 
bearing  from  S  . 

Optional:  Display  range  from  S, 


Repeat  Step  4  or  5  or  Step  5 
only  for  a  second  (subsequent) 
time,  then  proceed  to  Step  7. 

Compute:   Target  course 

Target  speed 

Display  Sensor   S   prompt. 

Target  bearing  from  S1   at 
latest  time. 

Target  range  from   S   at 
latest  time. 

Display  Sensor   S~   prompt. 

Target  bearing  from  S   at 
latest  time. 

Target  range  from  S„   at 
latest  time. 


course   ^ 
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speed   V 
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bearing  <J) 
range    p 
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9    (degrees) 
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fA 
fB 
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i 
c 
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D 
(R/S)' 
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(R/S) 

(R/S) 

(R/S) 
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none 
none 

course 
bearing 


R.  (n.mi.) 

l 
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\\>      (degrees) 
VT  (knots) 


Q        (degrees) 


R   (n.mi.) 


622  (degrees) 


r   (n.mi.) 


Step 

Instruction 

Input 

Key(s) 

Output 

8 

Compute  bearing  and  range  at  time 

t  (HH.MM) 

+ 

t   with  lead  distance  1    (n.mi.): 

I    (n.mi.) 

E 

a 

Display  Sensor   S..   prompt. 

1. 

b 

Target  bearing  from  S   at   t  . 

(R/S)* 

61£ 

(degrees) 

c 

Target  range  from  S   at   t  . 

(R/S) 

R£ 

(n.mi. ) 

d 

Display  Sensor   S   prompt. 

(R/S) 

2. 

e 

Target  bearing  from  S   at   t  . 

(R/S) 

62l 

(degrees) 

f 

Target  range  from  S.  at  t  . 

(R/S) 

Vl 

(n.mi.) 

9 

Repeat  from  Step  3  or  from 
Step  6  as  required. 

Note:   The  (R/S)  function  is  required  when  using  the  HP-67  mode.   This  output 
is  automatically  printed  on  the  HP-97. 


2 .   Sample  Problem 

a.   The  Primary  Sensor   S,   is  traveling  on  a  course  of  210°  at 
10  knots. 


b.  At  the  time  of  the  first  contact  sensor,  S_   is  115°  and 
3.5  n .mi .  from   S, . 

c.  At  1200  hours  the  first  contact  is  at  245°  from   S.   and 
260°  from   S_.   How  far  is  the  contact  from   S,   and   S?? 
(Ans . :  8  n.mi.  from   S,   and  10  n.mi.  from   S?.) 

d.  At  the  next  time  mark  sensor   S„   is  100°  and  5.0  n.mi.  from 

V 

e.  This  next  time  mark  is  at  1230  hours  with  the  contact  at  160' 
from   S,   and  239°  from   S?. 

f.  Estimate  the  course  and  speed  of  the  contact. 
(Ans.:   126°  and  14  knots.) 


g.   What  is  the  bearing  and  range  of  the  contact  from   S,   at 
1230  hours?   (Ans.:   160°  and  3  n.mi.) 
From   S2?     (Ans.:   239°  and  4  n.mi.) 


h.   Estimate  the  bearing  and  range  of  the  contact  from  S, 
Sp   at  1245  hours  with  a  lead  distance  of  3.5  n.mi. 
(S,   Ans.:   137°  and  10  n.mi.) 
(S?   Ans.:   164°  and   7  n.mi.) 


and 
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3 .   Program  Storage  Allocation  and  Program  Listing 
Registers : 


RO:   X2 

T 
Rl:   Y2 

R2:   t2 

R3:   At 


R4:   0 


R5:   6 


12 


22 


R6:   R. 


R7 


R8:   t, 


R9 :  I,    R  9      and   r 


SO: 

Xl 

SI: 

T 
Yl 

S2: 

4 

S3: 

S4: 

Zx 

S5: 

Used 

S6: 

£y 

S7: 

Used 

S8: 

Used 

S9: 

Used 

A: 


P 


i 
B:   (j,i 

C:   V 
s 

D:  \i> 
v  s 

E:   VT 
I:   4^ 


Initial  Flag  Status  and  Use: 
0:  OFF,  HP6  7  or  HP9  7  output 
1:   OFF,  Unused 


2:   OFF,  Used  for   t    option 
3:   OFF,  Unused 


Display  Status:   DSPO ,  FIX,  DEG 

User  Control  Keys : 
A:  ty      +  V 


B:   cb.  i  p. 

C:   t.  +  6. .  +  eo. 
l     li     2i 


a:  HP-67  output  mode 
b:  HP-97  output  mode 
c.   Unused 


D:   Compute  \p    ,  V   and  position   d:   Unused 


E:   t£  t  I 


e.   Unused 
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E.   Geometric  Analysis 

1 .   Static  Geometry 

Let   R.  =  (9,.,R.)   denote  the  bearing  and  range  of  the  target 
from  the  reference  (primary)  sensor   S,   at  time   t.,  and  let 
r.  =  (6   ,r.)   denote  the  bearing  and  range  of  the  target  from 
the  secondary  sensor   S„   at  time   t. ,  i  =  1,2,  where   t,  <  t„. 
Let   p.  =  (<}).,  p.)   denote  the  bearing  and  range  of   S9   from   S, 
at  time   t. .   The  static  geometry  for  some  fixed  time   t.   is 
depicted  in  Figure  1 . 

From  Figure  1  we  see  that 


R.  =  p.  +  r.  .  (1) 

li     l 


By  equating  the  rectangular  components  of  Equation  (1)  we  have 


R.  cos  9,  .  =  p  cos  d)  +  r .  cos  9_.               (2a) 

i       li  i       2i 

and 

R.  sin  9,  .  =  n  sin  $    +   r .  sin  90.  .             (2b) 

i       li    M  Y     l       2i 


Equations  (2)  are  two  equations  in  the  two  unknown  ranges   R. 
and   r..   Solving  this  system  of  equations  we  obtain 


sin(62i  -  4»i) 

R.  =  P.  — i — rz5 3 r        for  any   i,       (3) 

l     l  sin(^21  -  9.,)  xi 

and 
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N 


FIGURE  1.   The  Relative  Sensor  and  Target  Geometry 


at  Time   t . . 
1 
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sin(eli    -   <j>i) 
r.    =    p.    for   any      i    .  (4) 

1  sin(e2i  -  e  ) 


At  any  time   t.   the  target  range   R.   from  sensor   S,   and  the 

target  range   r.   from  sensor   S~  may  be  computed  from  Equations 

(3)  and  (4),  respectively.   Thus  R.   and   r.   are  determined  at 
any  time   t . . 


2 .   Dynamic  Geometry 

Let   V   =  (\\)    ,V  )   denote  the  course  and  speed  of  the 
primary  sensor   S, ,  and  let  V  =  (ty-,    V  )   denote  the  unknown 
course  and  speed  of  the  target.   Let   At  =  t~  -  t,  >  0   be  the 
time  between  first  and  second  observations  of  the  target.   The 
absolute  motion  of  sensors  and  the  target  is  depicted  in  Figure  2 
From  Figure  2  it  is  evident  that  one  of  the  many  vectorial 
relationships  is 


R,  +  Vm  At  =  V   At  +  R_  .  (5) 

IT        s        2 


The  target  course  and  speed  vector  V_   is  then  found  to  be 


^T  "  ^s  +  jt    («2  "  V  •  (6) 
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T(t  ) 


FIGURE  2.   Motion  of  Sensors   S,   and   S2   and  of  the  Target 
T   from  Time   t,   to  Time  t^. 
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3 .  Lead  Distance  Geometry 

If,  at  some  time   t   (t   >  t  ),  it  is  desired  to  lead 
the  target  on  its  track  by  a  distance  I,    then  the  bearing   0, p 
and  range   Rp   to  this  position  from  the  primary  sensor   S,   is 
obtained  by  converting  the  vector  [\p    ,  V  (t   -  t_)  +  I]       to 
rectangular  coordinates  and  adding  it  to  the  rectangular  form 
of  the  position  vector   R~   (see  Figure  3) .   The  resulting  vector 
is  then  converted  to  polar  coordinates  to  obtain  the  vector 
(6,0,  R0).   The  predicted  bearing  and  range   r    of  the  target 
from  the  secondary  sensor   S?   is  computed  from 


r£  =  R£  -  p   .  (7) 
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Position  of  target 
at  time   t„ 


Position  of  sensor 
S,  at  time  t 


Target  estimated 

position  at  time  t 

with  respect  to 

sensor 

at  time   t 


Si's  position 


Estimated  position  leading 
target  by  a  distance  £  at 
time   t^  with  respect  to 
the  position  of  sensor   Si 
at  time   t  . 


FIGURE  3.   Target  Lead  Distance  Geometry 
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